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Abstract. The stochastic systems without detailed balance are common in various 
chemical reaction systems, such as metabolic network systems. In studies of these 
systems, the concept of potential landscape is useful. However, what are the sufficient 
and necessary conditions of the existence of the potential function is still an open 
problem. Use Hodge decomposition theorem in differential form theory, we focus on the 
general chemical Langevin equations, which reflect complex chemical reaction systems. 
We analysis the conditions for the existence of potential landscape of the systems. By 
mapping the stochastic differential equations to a Hamiltonian mechanical system, we 
obtain the Fokker-Planck equation of the chemical reaction systems. The obtained 
Fokker-Planck equation can be used in further studies of other steady properties of 
complex chemical reaction systems, such as their steady state entropies. 
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The nonequilibrium behaviors of a great amount of chemical reaction networks can 
be described by stochastic differential equations. For example, the kinetics of a simple 
chemical chain reaction, 

^A^B^, (1) 

can be described by a set of stochastic differential equations, namely, chemical Langevin 
equations [HE} [3], 

dxi y/kixi 



dx 2 Vhxi . . 

— = k lXl - k 2 x 2 + — =-C 2 , (2) 



where Xi denotes the concentrate of ith species in the chemical reaction system of 
interest, here, x\ = [A] and x 2 = [B] ("[ ]" denotes concentrate of the specie). The 
fcj are rate constants for corresponding reaction steps, % = 0, 1, 2 stands for the reactions 
in Eq. (JTJ from left to right. The Q is the closed system's volume, Q, i — 1,2 refers to 
stochastic varible. 

We use the state variables qi = lnxj, instead of X{ G 1Z + to make sure g.; G 1Z. Here, 
we do not require the detailed balance condition to be satisfied in our complex chemical 
reaction system. In general, chemical Langevin equations are nonlinear stochastic 
differential equations. 

It is convenient to rewrite Eq. (J2]) to a concise form. From the point of view of 
tensor field theory (4], Eq. (|2j) can be written as, 

q a = w a (q) + h ab (q)( b , (3) 

here, x a , ( a and w a (q) are all (0, 1) tensor fields, or 1— form in TZ n , in the language 
of differential form theory. The h a b is a symmetric (0, 2) tensor field, h a t> = hb a . 
Here, we have used Penrose's abstract index notation [1], for example, ( a = (e M ) a , 
fj, = 1,2, ... ,n are the components of the ( a in the basis {(e M ) a }. In particular, if 
use the special basis {(e M ) a = (dx M ) a }, which are the dual basis of natural basis {(<9 M ) a }, 
Ca = Cm (.dx^) a . The following discussions can be applied to the systems with a large 
number of species, or mathematically, the systems have high-dimensional state spaces 
(manifolds) lZ n , with coordinates {q^ 1 }, \i = 1,2, . . .n. Note that Eqs. (j2J) have broad 
applications in the researches of physical, chemical, ecological, and financial systems, 
Eq. (j2J) is not merely restricted to chemical reaction systems [51 E]. 

In previous studies [8j [9], Eq. (J2]) was treated in the mathematical language of 
vectors and matrices. However, in order to effectively exploit the power results of 
modern mathematics, in present work, we will treat the chemical Langevin equations in 
a different way. In the studies on the existence of nonequilibrium processes described 
by stochastic different equations, Ao et al. [SJ [9], [TOl [11] introduced an "axillary matrix" 
M(g) = S(g) + A(g), and S(T)(g) is symmetric (antisymmetric) part of matrix M(g). 
The meaning of the the approach can be much clearer using the language of differential 
form theory. By introducing an axillary tensor field M a b , original differential 1-form 
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w a in Eq. [3]) is converted to another differential 1-form w a = M a b wv The existence of 
potential means 

«5a = (d/) a , feC°°(U), (4) 

and U is open in 1Z n . The necessary conditions of Eq. (J4]) can be expressed as following 
equations [TJ, 

g? = g, 0^=l,2,...,n). (5) 

The Eq. ([S]) restricts the choices of axillary tensor field M a b . It is worth noting that 
Eq. (J5J) can not ensure that the potential function is / G C°°(U). 

In present work, we directly decompose the differential 1-form w a using Hodge 
decomposition theorem, instead of introducing an axillary tensor field. The Hodge 
decomposition theorem is the extension of Helmholtz decomposition theorem in 
traditional vector analysis of TZ 3 . It tells us, for any differential p-form w p , (p < n) 
in lZ n , it can be uniquely decomposed into three differential forms, 

w p (q) = daip-^q) + 6/3 p+1 (q) + r/ p (q), (6) 

Here, the d and 5 denote external derivative and codifferential operators 
respectively [7J. The a is a differential (p — 1)— form, (3 is a differential (p + 1) — 
form, and 7 is a harmonic p— form, i.e., A7 = (dS + 5d) , -f(q) = 0, and "A" is the 
Laplace-Beltrami operator. 

In particular, Hodge decomposition of differential 1-form in 7Z 3 is equivalent to 
Helmholtz decomposition in traditional vector analysis: given a vector V, it can be 
uniquely decomposed asV A = V- + VxA, therefore Hodge decomposition can also 
refer to Hodge-Helmholtz decomposition. 

For our problem of interest, 1-form w a can be decomposed as w a = (d<p) a +(5^) a + r y a , 
with the being a scalar field (0-form), (3 a b being a 2-form (antisymmetric (0,2) tensor 
field) and 7 a being a harmonic 1-form, i.e., (d8 + 5 d)j a = 0. 

In present problem, the index of the tensor field can be lowered and raised by 
metric tensor in lZ n , g a b = S a b. Using Hodge decomposition theorem, after certain 
mathematical deriving, we obtain a stochastic differential equations with a form, 

ii^-^ + ^ + Mti + hM?, z = l,2,...,n. (7) 

Here, "<&" denotes time derivative of qi(t). The fij is antisymmetric, i.e., = —fji- 
The constraint on Ai, d 2 Ai(q) = implies Ai(q) is a harmonic function. We will 
use the properties of harmonic function later. Here, we have already adopted Einstein 
summation convention. In Ref. [8], with axillary matrix 5(q*) and A(qt), the so called 
^-decomposition need further constraints, and the constrained ^-decomposition is then 
called "gauged" <j)— decomposition. In our approach, the decomposition is definitely 
unique, which is guaranteed by Hodge decomposition theorem. 

To obtain the Fokker-Planck equation of our complex chemical reaction systems, 
a powerful method is the mapping from chemical Langevin equations to a Hamiltonian 
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mechanical system [SJ [HI [ED] • The deterministic forces and stochastic forces in Eq. (0) 
implies there are two different time scales in the system, which describe fast stochastic 
motions and slower chemical kinetic evolutions respectively. As discussed in Ref. [8], 
slow motion can be thought of as the consequence of the inertial of the chemical reaction 
systems, which can be described by an inertial mass m in the theory. 

We improved the approaches presented in Ref. [9], map the chemical Langevin 
equations to a Hamiltonian system which has a 2n— dimensional phase space, then obtain 
the Klein-Kramers equation for time evolution of probability distribution function. From 
now on, we will not use previous mathematical language of differential form, but directly 
use the components of differential forms, such as pi, fij, etc. To construct the 2n- 
dimensional phase space, or so called "enlarged state space", we introduce dynamic 



momentum pi, 



let 

4i =Pi/m, 



1,2, 



n. 



The other n Hamiltonian equations are constructed based on Eq. ((7j), 

dfij(q) 



Pi d(f)(q) 

Pi = TTT - 

m oq % 



dq j 



MQ) + h tj (q)C 3 



(8) 



(9) 



There is no Itd-Stratonovich dilemma in the connection between the stochastic 
differential equations, Eq. ([2]) and the Hamiltonian equations, Eq. fl8]) and Eq. (Q. 

We will truncate the results at the order of 0(l/m), finally integrate out the 
freedom of dynamic momentum p, and obtain the time evolution of p(q,t). For 
simplification, we assume a noise in Eq. ([9]) is a standard Gaussian white noise with m 
independent components Q, i = 1, 2, . . . , m. Define £j(g) = hij 
1,2, ... ,m, we have (&) = 0, and 

(&(%(r)> = 2D(q)6 ij 6(t-r 



1,2, 



,n, J 



i,j = 1,2, 



,71. 



(10) 



Here, the number of stochastic variables m can be different than that of state variables, 
the matrix D(q) is diffusion matrix. 

The generalized Fokker-Planck equation, so called Klein-Kramers equation can be 
obtained from Eq. (jHJ) and Eq. (jHJ), 

'Pi, 
_m 



dtpi(q,P,t) = dp 



+ d qi (j) (q) - d qk f lk - Ai{q) + Aj(g)<9j 



pj 



pi(q,p,t) 



Pi 



m 



which can be rewritten as, 



dt + V ~d q , 
m 



Gi(q)d Pi - d p 



Pi 



in 



Pi{q,p,t) = o. 



(12) 



Here, G^q) = d q J ik (q) + A^q) - (f) qi (q), and <p qz = d qi (j) (q). 

To get Fokker-Planck equation, we regroup the p— and q— derivatives in Eq. (fT2~ 
into two operators L\ and L 2 , 

P 



0„ 



D-d 



Pi » 

On 



m 



m 

i<t>M 



Mq) - d g k fik}di 



Pi ■ 



(13) 
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and Eq. (TT2"j) becomes, 

d t p(q,p,t) = {U. + L 2 ) Pl (q,p,t). (14) 

Then use Gardiner's standard projection operator method to eliminate the fast 
degrees of freedom of q implied in the zero mass limit, or "over-dumping" limit. We 
introduce a projection operator P, as 

Ph(q, p, t) = \ exp (-Pi^LPi) [ h (p\ q\ t)dV- (15) 

(V27rmdet Dj \ 2m J J 

here h(q,p,t) is an arbitrary function, and detD is the determinant of the diffusion 
matrix D(q). 

Projection operator is an idempotent operator, P 2 = P. Notice the fact, 
^exp 

my \ 2m 



D l0 d P] + ^ ) exp ( -^KtL ) = 0, (16) 



obviously there is an identity, L\P = 0. For an arbitrary ^-independent Cj(g), the 
following identity holds, 

w , )exp (_^)._I PA( , )exp (_^). (17) 

Thus, ip = Pi c i(<l) ex P (~ PkI 2m * s ^ ne eigenfunction of operator L\, with the 
eigenvalue being — 1/m, and the inverse operator L^ 1 satisfies, 

= -mV>. (18) 
Apply projection operator P to the phase space probability distribution function 

p(q,p,t) 

hi[q ' p ' t)= ^^m , " p (~ P ~ 1 ^ 1 ) p[q ' t] ' v(q ' p ' t) ' (19) 



with p(q,t) = J dp pi(q,p,t) being the state space probability distribution. 
It has been proved [5] 

d t v(q,p,t) = -PL 2 L^L 2 v + O (y/m) . (20) 

Substitute of explicit forms of operators L\ and L 2 to Eq. ([20]) . we have 

- PL 2 L^L 2 v = d gi [D a d gi + Mq) + ^0(g)] v(g,p,t)- (21) 

Finally, by integrating out the degrees of freedom of p in p(q,p,t), we get the 
multivariate Fokker-Planck equation with the form, 

d t p(q,t) = d qi [aa + Mq) + 9 q Mq)] pM- (22) 

We have known Ai(q) is a harmonic function. Assuming it is bounded at the state 
space TZ n , it must be constant, (denoted as Ai(q) = Oj), according to the Liouville's 
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theorem in harmonic function theory [T5]. With this assumption, our Fokker-Planck 
equation can be further reduced to 



Here, <f)(q, t) = t) + q\ is the potential landscape of the chemical reaction systems, 
therefore the "force" has a form of 



To conclude, we have used Hodge decomposition approach to study the potential 
landscape for chemical reaction systems described by chemical Langevin equations, 
which are stochastic differential equations. We do not require detailed balance condition 
to be held in the systems. Hodge decomposition theorem is very convenient to deal with 
the problem related to high-dimensional state space. We map the chemical Langevin 
equations to a Hamiltonian mechanical system, and find the Fokker-Planck equation 
for chemical Langevin equations. With certain assumptions, there does exist potential 
landscape, the potential can also be time-dependent. Our work can be used in the 
studies of other steady state properties of the chemical reaction systems, such as their 
steady state entropies. 
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d t p(q,t) = d qi Dij(q)d qj +d qi (f) p(q,t). 



(23) 



fi = -[Mq) + d«<KQ)] 



[<k + 9 qi 4>(q)} 



d qi (f>(q)- 



(24) 
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